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^ ' Abstract 
> ■ 

, In the preceding paper arXiv : 0802.3252 [quant-ph] we treated a model given by a 

' master equation with generahzed Lindblad form, and examined the algebraic structure 

■ 

■ related to some Lie algebras and constructed an approximate solution. 

oo ; 

, In this paper we apply a unitary transformation by the squeezing operator to the 

. ' master equation. Then the generalized Lindblad form is tranformed to the usual Lindblad 

X ■ 

\^ • one, while the (original) Hamiltonian is tranformed to somewhat complicated one. 

As a result we have two different representations based on the Lie algebra sn(l,l). 
We examine new algebraic structure and construct some approximate solution. 

Quantum Computation (Computer) is one of main subjects in Quantum Physics. To realize 
it we must overcome severe problems arising from Decoherence, so we need to study Quantum 
Open System to control decoherence (if possible). 
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This paper is a series of [T], [2] and [3] and we continue to study dynamics of a quantum 
open system. First we explain our purpose in a short manner. See |1] as a general introduction 
to this subject. 

We consider a quantum open system S coupled to the environment E. Then the total 
system S + E is described by the Hamiltonian 

Hs+E = Hs®1e + 1s® He + Hi 

where Hs, He are respectively the Hamiltonians of the system and environment, and Hj is the 
Hamiltonian of the interaction. 

Then under several assumptions (see [1]) the reduced dynamics of the system (which is not 
unitary !) is given by the Master Equation 

^^p=-^[Hs,p]-Vip) (1) 

with the dissipator being the usual Lindblad form 

V{p) = ^Y. {"^^^^P + A- - 2^.P^]) • (2) 

{j} 

Here p = p{t) is the density operator (or matrix) of the system. 

Similarly, the master equation of quantum damped harmonic oscillator (see [1], Section 3.4.6 
or also [3]) is given by 

—p = —i[ua^a, p] — ^ {a) ap + po) a — 2apa^) — — {aa) p + paa) — 2a) pa) , (3) 

where a and are the annihilation and creation operators of the system (for example, an 
electro-magnetic field mode in a cavity), and /i, p [p > v > Q) are some real constants 
depending on the system (for example, a damping rate of the cavity mode). 

In p] the general solution was given in the operator algebra level. This is a very 
important step. Moreover, in [T] we treated the master equation with generalized Lindblad (or 
Kossakowski-Lindblad [U]) form given by 

—p = —i[ua^a, p] — ^ [a^ap + pa^a — 2apa^) — — (aa^p + paa^ — 2a^ pa) 

(aV + pa^ - 2apa) - ^ ((a^)V + p(a^)^ - 2a^pa^) (4) 
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where k is a complex constant satisfying the condition fiu > which ensures the positivity. 
Then we examined the algebraic structure related to the Lie algebras su{l, 1) and su(2), and 
constructed some approximate solutions by use of it. 

In this paper we apply a unitary transformation by the squeezing operator to (jlj), which 
was developed by An et al. [Z] , [H] □• Then the generalized Lindblad form is tranformed to the 
usual Lindblad one, while the (original) Hamiltonian is tranformed to somewhat complicated 
one. We examine new algebraic structure in detail. 

The squeezing operator S{e) is given by 



(5) 



^(e) = exp Q (e(at)2 - ea^)^ , e = \e\e'^ G C 

(|e| and (p are chosen later on) and it has well-known remarkable property 

S{e)aS{ey = c{\e\)a - e'^s(|e|)a^ S{e)a^S{ey = c{\e\)a^ - e-'^s(|e|)a (6) 

where c(|e|) = cosh(|e|) and s(|e|) = sinh(|e|) for simphcity, see for example [9]. 
We act the squeezing operator to the master equation (jH) as adjoint action 

^Ps = -i[^as^as,ps] 

{as^a-sPs + PsO'S^'as - 2aspsas^) - ^ (a^asVs + Psasas^ - 2as^ psas) 
"2 [^s^Ps + Psas^ - 2aspsas) - - ((as"^)^^ + Psias^'f - 2as^ psas^) (7) 

where ps = SpS'' and as = SaS'' . Now let us calculate the right hand side. 
It is not difficult to check 

As = as^asps + Ps^s^as - 2aspsas^ 

= (? [^a)aps + pso)a — 2apsa)^ + [^aa) ps + psao) — 2a) pso] 

-e-'^cs {a^ps + psa^ - 2apsa] - e"^cs {{a^fps + Psia'^f - 2a^ psa^] (8) 



^ In the papers the positivity condition > |kP has been neglected 



3 



and 



Bs = agas^'ps + Psaso-s'' - '^ag^ psag 

= [^a)aps + pso'a — 2apsa'^ + {aa^ ps + psaa' — 2a^ psa^ 

-e-'^cs {a^ps + psa" - 2apsa} - e'^cs {{a^)^ps + psia^f - 2a^psa^} 



and 



Cs = alps + psal - 2aspsas 

= —e^'^cs {a'aps + psa'a — 2apsa'} — e^'^cs {aa' ps + psaa) — 2a) pso^ 
{a^ps + Pso? - 2apsa] + e^'h'' {{o)fps + Ps{a)f - 2o) Pso^] 



and 



Ds = {as^fps + Ps{as^f - 2as^psas^ 

= —e~^'^cs \^a^aps + Ps^^a — 2apsa^} — e~^'^cs \^aa^ ps + psaa) — SaVs^} 
+e-2^^s2 {a^ps + psa^ - 2apsa} + {{a^fps + P5(a^)' - 20^5^^} 

where c = c(|e|) and s = s(|e|) for simplicity. 

Therefore by (jH]) ~ ( fTTl) the generahzed Lindblad form in ([7]) becomes 

~ {-pe-'^cs - ue-'^cs + kc^ + e-^''''Rs^) C-^ {-pe'^cs - vd'^cs + Rc^ + e^^-^Ks^ 



where A = As=i et al. 

Now we can remove the terms ffT^ if we choose e suitably. In fact, we set 

K = \K\e-'^ = ke~'^ 

then we have a quadratic equation 

- (/i + v)cs + k{c^ + s^) = 4==^ - (/X + + A; = 



for t = ^ = tanh(|e|). The discriminant V is 

V = {^ + uY - 4:k^ = (/U - + 4(/iz/ - A;2) > 

by the positivity condition fiu > k"^, so there is a solution on t on |e|). 

A comment is in order. We cannot remove the second term in (fT2l) against the result in 
jH], [Zj by the positivity condition. We leave it to readers. 

Therefore, the generalized Lindblad form becomes 

/ic^ + z/s^ — 2kcs . I X /is^ + uc^ — 2kcs , ^ x x . 
^ \a}aps + ps<ra- 2apsa^) [aa^ Ps + Psaa^ - Psa) ■ 

(15) 

Next let us calculate the Hamiltonian in ([7]), which is easy to see 

- i[ujas^as,ps] = {(c^ + s'^)a^a - e"''cs{a^f - e-'^csa^} , ps]. (16) 
As a result the master equation (IT!) becomes 

= {{c^ + s^)a^a - e"^cs{a^y - e-''^csa^} , ps] 



/ic^ + us^ — 2kcs 
2 

/is^ + uc^ — 2kcs 



{a)aps + psa)a — 2apsa)) 

{aa) Ps + psCLo) — 2a) ps<x) (17) 



under the squeezing operator S{e) with some e = lele*"^. 

Next we apply the method in [1], [2] to the equation (fTTI) to make the algebraic structure 
clearer. Before that let us make some mathematical preparations. 

A matrix representation of a and on the usual Fock space 

T = Vectc{|0), |1), |2), |3), ■ ■ • }; \n) = 

is given by 
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a = e 



id 



1 


73 


/n 



at = e"*^ 





1 
V2 





\ 



J 



where e*^ is some phase. Note that ao) = a)a + l = N+l. 
For a matrix X = (xij) G M{T) 



^ Xii Xi2 ■ ■ ■ ^ 



X 



X21 X22 X23 ■ ■ ■ 
X31 X32 X33 ■ ■ ■ 

V • • • ■• 



we correspond to the vector X E J^^^^c^ as 



(18) 



(19) 



^ — i^ij) ^ — (3^115 2^12) 3:13, 2:21, 2:22, a;23! ^31) ^32; 2:33, 



r (20) 



where T means the transpose. The following formula 



AXB = {A®B' )X 



(21) 



holds for A,B,X e M{T). 
If we set 

K^ = ]^{N ®1-1®N), 



(22) 
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where A^"^ = A^, then we have 

[i^3, I<+] = I<+, [K3, K^] = -K^, K^] = -2Ks. (23) 

If we also set 

K3 = ^{N®1 + 1^N + 1^1), k+ = a)®a^, k_ = a®{a)f, (24) 

we have 

[^3, k+] = [ks, k_] = -K., [k+, k^] = -2ks. (25) 

Namely, {K^, K^, K_} and {K^, K_^., K_} are a set of generators of 5-^(1, 1) algebra. 
Now the equation (ITTj) can be written as 

= J -2iu{c'^ + s'^)K3 + 2iue'^csK+ + 2iue-''^csK_ + ^^-^1 ® 1 

- ((/i + z/)(c^ + s^) - Akcs) k-s + {fis'^ + uc^ - 2kcs)k+ + {fic^ + us^ - 2kcs)k_^ ps 
^'""l^l + A + Mps] (26) 



2 

A = -2iu{c^ + s^)K3 + 2iuje''''csK+ + 2iuje-''^csK_, (27) 
A = -{{p + v){c^ + s^) - Akcs) ks + {ps^ + vc^ - 2kcs)k+ + (/ic^ + vs^ - 2kcs)k^ (28) 

by use of (1^ and generators fl2^ . fl^ . The solution is formally given by 

ps{t) = e^*e*(-^+-^)ps(0). (29) 

Now let us calculate the commutators of {A'3, K^, K_} and {K3, K^, K_} : 

[K3,i^+] = [K3,k^] = [K,,k,]=0 

and 

[K+, k^] = -at ® {a^, [K^, K_] = -at ® (at)^, + , 7^-3] = -^{(«^)' ® 1 + 1 ® ((at)')^}, 
=a®a^, = a®a^, [K^, k^] = ^{a^ ® 1 + 1 ® {aY}. 

Since {i^3, -f^-} and {K^, don't commute it is reasonable to assume 

ps{t) ^ e'^W-^psiO) (30) 



as the first approximation. 

The disentanghng formula (which is well-known) for the system {L3, L+, L„} based on the 
Lie algebra su{l,l) is given by 



with 



^;^sinh {Wa^ - be) 



^ ^ cosh (t^/a'^ — bcj — sinh (ty/ a? — bcj ' 

F{t) = cosh (tV — be] -j^=^= sinh ( a? — he] , 

^ ^ \/a? — be ^ ^ 

, ,f , sinh f — 6c) 

S(t) = . , -^^X ^ / , (32) 

cosh [ty — be) — sinh [ty a? — he) 

See for example [9], p^ . 

From this formula we can calculate e*"^ and e*-^ in (l30l) as follows. 

gM ^ ^G{t)K+^-2\og{F{t))Ks^E{t)K- ^^^^ 



with 



G'(t) 



^;^sinh (tVa2 - 6c) 



and 



cosh {t^/ — be) — sinh (^t^/ a? — he) ' 

F(t) = cosh (tV — bc\ , sinh {t\J a? — he 

sinh (tV — 6c) 
cosh — 6c) — ^^2_i,f, sinh — 6c) 



a = -ia;(c2 + s2), 6 = 2zcue^'^cs, c = 2zcje-^'^cs. (35) 



Then it is easy to see \J — he = iu and 

^/ N 2ze*'^cs sin(ta;) ^, , , , ., , on . / n 

G(t) = - ^/-J. V ^(^) = cos(t^) + z(c2 + s2 sm(to; , 

cos(rc<jj + «(c^ + s^j sm(ia;) 

^ = 2^e-''^cgsin(ta;) 

^ ' cos(tcj) + 2(c2 + s2) sin(tcj) ■ ^ ' 

Next, 
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with 



and 



^ ^ cosh (tVtt^ — be) — sinh (t\^ a? — be) ' 

F(t) = cosh (tV — he\ , sinh ( — he] , 

sinh {t\/ a? — be) 
cosh (t^/a'^ — be) — sinh (t V — be) 



(/i + ^)(c^ + g^) 4A;cg ^ 6 = Ais^ + j.c2 - 2A;c5, c = /xc^ + i^s^ - 2te. (39) 



\/ — be = — . 
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Then it is not difficuh to see 



However, by use of it we cannot simphfy fl38l) hke 
Therefore our approximate solution (l30l) becomes 



ps{t) ^ e^*e^W^+e-2i°s(^W)^3ei^Wi^-eGWi?+e-2iog(F(t))X3e^WX-^^(0) (40) 



and we restore this form to the usual one by use of (12T1) . The result is 
ps{t) ^e^*exp (^^{a^y^ {exp (- log(F(t))iV) x 

^exp (^^a^) 0(t)exp (-^a^) } exp (log(F(t))iV)| exp (^-^{a^^ (41) 
and 



E ^«>^(0)(«^)™| (-log(F(t))iv)} a". (42) 



m! 



This is indeed complicated. Compare this with the corresponding one in p. 

In this paper we revisited the quantum damped harmonic oscillator with generalized Lind- 
blad form and applied the unitary transformation by the squeezing operator to the master 
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equation, and examined the new algebraic structure and next constructed some approximate 
solution in the operator algebra level. 

The model is very important to understand several phenomena related to quantum open 
systems, so the general solution is indeed required. To obtain it (like in [2]) is almost impossible 
at the present time. 

Lastly, we conclude the paper by stating our motivation. We are studying a model of quan- 
tum computation (computer) based on Cavity QED (see [H] and [12]), so in order to construct 
a more realistic model of (robust) quantum computer we have to study severe problems coming 
from decoherence. 

For example, we have to study the quantum damped Jaynes-Cummings model (in our 
terminology) whose phenomenological master equation for the density operator is given by 

—p = —i[Hjc, p] — — {a)ap + pa) a — 2apa^) (aa^ + paa^ — 2a^ pa^ , (43) 

where Hjc is the well-known Jaynes-Cummings Hamiltonian given by 

Hjc = yfTs ® 1 + uoh ® a"^a + (a+ ® a + a_ ® a^) (44) 





with 

0"+ = ' ^ = ' = , I2 = 

yooy \^ I \^ ^ 

Note that p G M(2; C) ® M{T) = M(2; M(J^)), where M{T) is the set of all operators on the 
Fock space JF. See for example [13], [H]. 

Furthermore, it may be possible to treat the generalized master equation given by 

—p = —i[Hjc, p] — ^ [a^ ap + pa'^a — 2apa^^ — — (^aa^ p + paa^ — 2a^ pa^ 

— ^ (aV + pa^ - '2apa) - ^ ((a"^)V + pi^^f - '2a^ pa^) (45) 

with the condition pu > |Kp similarly in this paper. 

These equations ( (I43p . 0451) ) are very hard to solve in the operator algebra level, so even 
constructing approximate solutions is not easy. This is our future task. 
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